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Summary. We briefly review the main aspects of (Poincare-Dulac) normal 
forms; we have a look at the non-uniqueness problem, and discuss one of the 
proposed ways to "further reduce" the normal forms. We also mention some 
convergence results. 

1 Poincare normal forms and their use 

Normal forms were introduced by Poincare as a tool to integrate nonlinear 
systems; by now we know this is in general impossible, but it turned out that 
the usefulness of normal forms goes well beyond integrability. 

An introduction to normal forms is provided e.g. in || |5 , 26 , 4l|, |63| ; 



see also (H], |33| . An introduction to normal forms for Hamiltonian systems is 
given in appendix 7 of Q. 

Their use is, beyond integrability issues, of three kinds: 

• Classification of dynamical systems; 

• Qualitative study of dynamical systems; 

• Quantitative study of dynamical systems. 

In this note I want to discuss the Poincare-Dulac approach, and an exten- 
sion thereof, focusing on the simpler setting, i.e. normal forms for a smooth 
dynamical system (equivalently, an autonomous first order ODE, a vector field) 
in R™, around an equilibrium point. For more general settings, see e.g. || S]. 

I consider a first order ODE (dynamical system) in R n having a fixed point 
in the origin, and written as a series of homogeneous polynomial terms: 



Ax + J2f k (x) (1) 



fc=i 
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where fk(ax) = a k+1 fk{x) (we say also G Vk, and write V — ©kL Vfc). We 
can of course think this is obtained from a general smooth^ equation x = f(x), 
with /(0) = 0, expanding itQ around x = 0. Actually, for most of our discussion, 
we could as well allow the series on the r.h.s. of (1) to be just a formal series 
(i.e. not require convergence in any neighbourhood of the origin). 

I will also assume that A is semisimpl^ for a discussion of the general case, 
see ^7|, |6J| and || [l9| Recall that if A is semisimple, it can be brought 
to diagonal form, and thus to be normal, [A, A + ] = (this will be the relevant 
property); it will thus be more economic to perform all computations in a basis 
where A is diagonal. 

The idea of Poincare was to eliminate nonlinear terms by a sequence of 
near-identity changes of coordinates, 

x = y + h m (y) (2) 

with h m G V m , and m taking successively the values m = 1, 2, .... 

With this, and writing D for the matrix of partial derivatives of h m , eq.(l) 
reads [I + D]y = f{y + h m (y)), which means 

oo oo 

y = [1 + D}- 1 My + h m {y)) := £ J k {y) (3) 

fc=0 fc=0 

By a simple computation, we obtain that fk(x) = fk(x) for k < m, and 
(see below for the notation) 

f m {x) = f m (x) - C (h m ) ; (4) 

the effect on higher order terms is more involved, and is usually not considered 
in detail, as higher order terms can be normalized at later stages. 

In eq.(4) above, Cq is a linear operator, defined in terms of the linear part 
A of / as 

£ o (V0 := (Ax-V)i> - Aip ; (5) 

notice that Cq : Vk — * Vk- Actually it is more convenient to describe this is 
terms of a general bracket 

:= {<p-V)il> - {^-V)ip (6) 

(this is just the representation of the commutator of the vector fields ip l di and 
ip l di in terms of their components); with this, Cq is just defined as Co{ip) := 
{Ax,tp}. 

1 By this we will always mean C°°. 

2 In this way we lose all the information "beyond all orders" in perturbation theory, such 
1/2 

as e 1 ' x terms; for this, see the works by Ecalle. Notice that if we deal with pcrturbativc 
power series expansions of smooth dynamical systems (vector fields), it would be natural in 
some aspect to consider C°° rather than analytic conjugacy. 

3 Let me stress, however, that a good part of the discussion conducted here will also apply 
to the case where A is not semisimple (at least provided it has nonzero semisimple part). 
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In this way we can eliminate all the terms which are in the range of Cq 
(note that if A is normal, Ker(£o) is complementary to Ran(£o) in V, and thus 
the whole normal form can be taken to be in Ker(£o); this means that the 
whole normal form commutes with the linear part of the system) . We can also 
introduce a scalar product^ in V, so that Ran(£ ) = Ker^J). 

We can now also easily characterize (with A taken to diagonal form) the 
terms in Ker(£j); assume Ai, A„ are the eigenvalues of A, and vi, v„ the 
corresponding eigenvectors (the computations are immediate if we take a basis 
of eigenvectors in R n ). Then a vector w = x M v r := x^ 1 ...x£ n v r is in Ker(£j) 
if and only if there are nonnegative integers (with \fj,\ := X^=i > -O 

such that 

n 

2^/x s A s = A r (7) 

s=l 

and Ker(£o") is spanned by such "resonant monomials". The equation (7) is 
called a resonance relation of order \fi\. 

Thus, following the Poincare procedure, we can always eliminate term by 
term all the nonlinear terms up to any desired order N, provided the A^ are 
nonresonant; note that this is generically the case. 

In the case where the A^ satisfy some resonance relation, i.e. when (7) 
admits some solution, we can always eliminate nonresonant terms, and arrive 
at a form of /, the Poincare- Dulac normal form, in which only resonant terms 
are present. That is, we arrive at coordinates in which (1) is written as[] 

oo 

x = Ax + J2dk(x) ; g k e Ker(£+) . (8) 

k=i 

It should be stressed that this result also holds for non-semisimple A [§3| . 

In this way we obtain a description of nonlinear systems (in normal form) 
"compatible" with a given linear part. 

In practice, the normalization will be carried out only up to some finite 
order N; one will then use perturbative arguments to infer that the dynamics 
of the full system is only a small modification - at least in a sufficiently small 
neighbourhood of the origin (and/or for sufficiently short times) - of that de- 
scribed by the truncated normal form. In this way, if we were able to study 
the most general (truncated) system in the form (8) for given A, we would in 
principles obtain a classification of the qualitative local behaviours of nonlinear 
systems around a fixed point with prescribed linearization A. 

Let me point out, after mentioning these "qualitative" uses of normal forms, 
the basis for its use in quantitative studies. Consider again (1), and denote by 
ji the first k with fk^O (so generically ju, = 1); rescale x as x — ex, divide by 

4 With the notation introduced below, the scalar product in V can be chosen to be the 
Bargmann scalar product ]47[ , (x^v r ,x v v s ) := S r ,s&^,v (mO- With this, Cq (tj}) = {A + x,ip}. 

5 More precisely, we can arrive to such a form up to order N, for any desired N. That is, 
we should add a remainder 7£jv °f arbitrarily high order N + 1 in this equation. 
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e, and drop the tilde for ease of notation, so that we have a system of the form 



x = Ax + er| M g(cc; e) ; (9) 

let x(t) — (fr e (xo,t) be its solution with initial datum x(0) = xq. If we approxi- 
mate this with the solution £(t) = &o(xo,t) to the system x = Ax obtained for 
e = 0, we are making an error 

Pe (t) ■= \x(t)-£(t)\ = |$ e (x ,t) - $o(*o,*)| (10) 

Assume that A has norm C, and that in a region B we have \g(x; e) | < M(e) 
(all of our estimates will be only valid when x and £ are both in B, i.e. for a 
time t depending on the exit time from B) . Then one obtains easilj]^] that 

Pe(t) < \e\^[e ct -l] . (11) 
This also means that p £ (t) < S for a time not smaller than 



to = h ln 



5C 

1 



(12) 



When we perform a normalization, with no resonance of order smaller than 
p, we are pushing nonlinearities to higher orders and thus increasing p in the 
previous estimates; that is, we are improving the accuracy of the linear esti- 
mate, or extending the time interval over which we have an error of the same 
"acceptable" size 6. 

It should be stressed that in changing coordinates we are also changing 
the value of the constant M; if the "small denominators" appearing due to the 
change of coordinates are too small, their effect will more than counterbalance 
the gain we obtain by pushing up the nonlinearities, and disrupt the perturbative 
expansion. 



2 The problem of convergence 

It should be noted that the discussion conducted above is purely formal: the 
changes of coordinates are constructed as a series which could - and in general 
will - be only formal. 

Here we want to consider the case of analytic systems, and for these we 
should study the convergence of the series defining the transformations. 

This is in general a hard problem (let us stress that it is difficult to prove 
convergence, but it is also difficult to prove divergence; see [5U5Q, and there 
are very few general results. Let me recall the main oncs[] (|[ 2C|27], |65|. Here 

6 Indeed dp E /dt < \A{x - f)| + \e^g(x, e)\ < Cp E + \e\^M (e). Thus the growth of p s (t) > 
is slower than that of a function r(t) obeying dr/dt = ar + b with r(0) = and a = C, 
b = \e\^M(e). The solution to this equation is r(t) = (b / a)[exp(at) — 1], 

7 It is remarkable that when we can guarantee convergence of the normalizing transforma- 
tion on the basis of the first four criteria below, we have an integrable normal form. 
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the function / : R" — > R™ denning the dynamical system x — f(x) should be 
thought of as / : C" -> C". 

• Poincare criterion: if the convex hull of the A, € C does not include zero, 
then the normalizing series converges in an open neighbourhood of the origin. 

• Siegel criterion: call the vector A = (Ax, A„) of type {C,v) if for each 
r = 1, ...,n and for any > 2 it satisfies | ^ s (/i s A s ) — A r | > C/\n\" '. If / is 
analytic and A is of type (C, v) for some positive numbers C and v, then / is 
analytically equivalent to its linear part. (Note that the set of vectors A which 
are not of type (C, v) for any C > has measure zero for v > (n — 2)/2; in the 
real case this holds for v > n — 1). 

• Pliss criterion: Pliss realized that Siegel condition can be weakened: call 
the vector A = (Ax, A„) of type P(K, v) if for each r = 1, n and for any 
Ia*I — K> it satisfies either the resonance condition 2 s (^i s A s ) — A r = 0, or 

^ s (/j s A s ) — A r | > I/ImT- If / is analytic, A is of type P{K,v) for some 
positive numbers K and v, and / is formally conjugated to its linear part, then 
/ is analytically conjugated to its linear part. 

• Bruno criterion: say that A satisfies "conditions" if J2T=i hi(l/wfc) < 
oo, where uik is the minimum of | ^2 S fi s X s — A r | on all the r = 1, n and all 
the < 2 fe not satisfying a resonance relation (condition co is weaker than 
the Siegel or Pliss condition ]l9|]). Say that / having linear part Ax satisfies 
"condition A" il^ f(x) = [1 + a(x)]Ax. Then, if A satisfies condition u and 
/ is formally conjugated to a normal form / satisfying condition A, then / is 
conjugated to / by an analytic transformation. (Condition A can be weakened 
for suitable classes of /'s). 

• BMW-C theory: assuming condition w for A, condition A can be replaced 
by suitable symmetry conditions on /. This theory was created by Markhashov 
[ [52"| and by Bruno and Walcher plj ; it was later extended by Cicogna |2f| and 
again by Walcher |39|. This matter is discussed in [^7j and in the SPT2001 
tutorial by Cicogna and Walcher p§] . See also the recent work by Stolovich |^9| 
which, together with many other things, contains results in this direction. 

• Normalization of an algebra of vector fields: one can consider the problem 
of simultaneously normalizing an algebra of (analytic) vector fields Xi around 
a common singular point. Vey |64[ ] proved that if they span an abelian algebra 
and their linearizations at the singular point are linearly independent, then they 
are simultaneously and analytically normalizable. Results in this direction are 
also contained in 0, E| |6Cj] ; see also Q . 

• Hamiltonian systems: in the case of Hamiltonian systems, we can deal 
directly with the Hamiltonian (i.e. a scalar function) rather than with the 
hamiltonian vector field this generates. Many of the results quoted above can 
be recast in this setting, but there have been investigations of the convergence 
in the hamiltonian case per se. In this context, we refer to the work of Russ- 
mann, Bruno, Vey, and Ito |l9|, 48, [54|, as well as to that of Giorgilli and 



collaborators |38| (also in connection with the study of effective stability and 

8 In geometrical terms, if / is orbitally equivalent to its linear part ^5tJ. This means the 
trajectories of / are the same as those of Ax. 
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KAM theory [JM E0| ) . I will not discuss these results here; see however the 



paper by Cicogna in || and |26|, |27| for brief discussions. 

• BCS theory: convergence of the normalizing transformation implies that 
the original vector field and its normal form are analytically conjugated. In 
many cases one would be satisfied with C°° conjugation [l3| (or even with C k 
conjugation: e.g., this is the only one making sense for (7*^ vector fields. For C k 
normal forms theory, see E^]). It was proved by Chen and by Sternberg |2^, Q 
that if the linear part A of the system is hyperbolic, then formal conjugacy 
implies C°° conjugacy (see also Jl^, |o[ Q). This result has been recently 
extended to equivariant (or reversible) systems, and also to the Hamiltonian 
framework, by Belitskii and Kopanskii E3|. 

• Orbital conjugacy: the reader should be warned that in many cases one 
is satisfied with orbital equivalence (two systems are orbitally equivalent if they 
have the same integral curves, independently of the time law they follow on 
these). This question is considered e.g. in || |]| |67j . 

In the following, I will come back to considering only formal transforma- 
tions; for a discussion of convergence see e.g. [^6[ |28|. 



3 Some remarks on Poincare normalization 

Now that I described the why and how of Poincare normalization, let me remark 
some simple points which are relevant for the following. 

• The change of coordinates needed to normalize / at order k is generated by 
hk solving Co(hk) — ito(fk), where ttq is the projection from V to the range 
of Co, and fk is the term obtained after the first (k — 1) normalizations. 
The solution to this "homological equation" for hk is unique up to a term 
in Ker(£o)- That is, if Shk = hk — hk G Ker(£o), then hk and hk define 
the same transformation in Vfc. 

• On the other hand, functions hk and hk as above have in general a dif- 
ferent effect on the f p with p > k. Thus the freedom in the choice of hk 
(for resonant systems) means that we can generate different normal forms 
starting from the same nonlinear system. 

• This shows at once that normal forms are somehow ambiguous, i.e. the 
normal forms classification is in general redundant: different normal forms 
(with the same linear part) can be actually conjugated. 

• It is maybe worth remarking explicitely that the higher order terms gen- 
erated by the change of coordinates x = y + hk(y) will in general also have 
a resonant part; this cannot be eliminated by the Poincare procedure at 
later stages^ 

9 Thus, in particular for resonant systems, it is possibly misleading to say that normalization 
is a way to "simplify" (a rather vague word) nonlinear systems: a system with very few terms 
can have infinitely many terms once taken into normal form. 



6 



It is clear that from the point of view of qualitative study and classification 
we would like to eliminate or at least reduce the redundancy in the normal 
forms classification. This would also have an obvious computational advantage: 
among normal forms corresponding to / we could choose the most convenient 
for our needs. 

Thus the problem of "further normalization" has been considered by several 
authorsR I will quote first of all Takens |6lJ], and Broer who in his thesis 
(advised by Takens) set the problem - and the theory of normal forms - in 
terms of filtrations of Lie algebras |]l6|, |l7| , |is}| . In this respect, one should also 
mention the early work of Kummer |5l[| . 

Broer's Lie algebra approach was also followed by van der Meer |53 and by 
Baider, Churchill and Sanders || JMLi * n particular they obtained a formulation 
of "unique normal forms" , see also |2J, [k| . Unfortunately, their results are very 
difficult to implement on concrete systems, even the simpler ones, and seems to 
be not applicable in concrete terms (on the other hand, Broer's approach is a 
guiding light for many questions in normal forms theory). 

Other, and more recent, attempts in this direction (some of them with 
different tools) include works by Algaba, Freire and Gamero, by Chen and Delia 
Dora, by Kokubu, Oka and Wang, by Ushiki, and by myself @, |||o[ H H 



Unfortunately, no systematic exposition and comparison of these different 
approaches is available; needless to say, this is not the place to attempt such a 
work, but maybe some reader of the present paper will endeavour a comprehen- 
sive discussion of this matter, which would surely be useful. 

The same problem of further normalization was also tackled, of course, by 
russian authors. One should quote in this respect at least the works by Belitskii 
[ [l2| and by Bruno pl|; these are also discussed and compared in the recent 
book |2(|. For russian work in the area see also, of course, ||. 

In the next sections I will present a simple wayPI to (in general, only par- 
tially) remove the degeneration of normal forms. This is based on ideas which 
I believe were first introduced by Broer (and recalled in a somewhere deformed 
way in the final section 6); the approach discussed below in sections 4 and 5 aims 
at a smaller generality than Broer's one, but is easily implemented in practice 
©[I- 

I would like to thank prof. Belitskii for pointing out, after a preliminary 
version of the present paper was circulated, that this approach is very much 
similar to the one proposed in his work [O . 



10 Actually already Dulac (see page 335 of |?^]) remarked that possibly a choice of the terms 
in Ker(£o) different from the one he was employing - i.e. Sh^ = - could have produced 

a simpler normal form. 

11 The reader should be warned that the discussion of it given in pp) uses different basic 
definitions; sec 136 1 
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4 Renormalization (further normalization) 



The approach I will discuss in this section (based on p5[ ; see also |26j| , and 
[ |36| for more recent discussions) is completely algorithmic, and not any more 
difficult to implement (by hand or by computer) than the standard Poincare one: 
actually it is just based on a careful exploitation of the freedom left by Poincare 
algorithm. As mentioned above, this is quite similar to the one proposed in |p"lf . 

The main idea is to control the effects of the change of variables generated 
by hk on higher order terms / p , p > k; as already mentioned, this change of 
variables will in general produce resonant terms. This should not be seen as a 
drawback, but rather as an advantage! By a judicous choice of hk, we can use 
this feature to eliminate higher order resonant terms. This "judicious choice" 
amounts to choosing the h^s as solutions to "higher homological equations". 

It is actually convenient to consider an improved version of Poincare trans- 
formations, i.e. Lie-Poincare ones; for some of the advantages of using these see 
[ [l5| , see also |66| . The idea is to define a vector field 

H ( 13 ) 

in R" and consider the flow <&(er;xo) generated by this. The change of coor- 
dinates is then given by x — $(l;y), which coincides with the standard one at 
lowest nontrivial order. 

The change of coordinates with generator h m maps (1) into a new system 
V = f(y) =Ay + Y,kLi fk(y), where f k (x) = f k (x) for k < m, while for k > m 
we have, with [a] the integer part of a, 

[fe/m] 

= £ ^rn(fk-s m ) ; (W) 
s=0 

here the operator H m is defined as H m (a) := {h m , a}, and means applying 
TL m for p times. Note that the expression (14) is just the translation to the 
present notation of the familiar Baker-Campbell-Hausdorff formula: indeed it 
states that the vector field X — pdi is changed into a vector field X = fdi 
given by X := [e aH Xe- aH ] a=r 

Having determined the effect of the Lie-Poincare transformation to all or- 
ders, we can try to control it. 

Let us consider again our system (1), i.e. x = f(x) = Ax + J2tLi fk(x)- 
We can normalize the quadratic term /i by solving the homological equation 
C {hi) = tt /i, and we ge10 

=fi-C (h 1 ) . (15) 

The higher order terms fk (k > 2) will of course also change with this first 
change of variables, but for the sake of notation we will not introduce indices 
to stress this fact, and call the new expressions of these again /2, /3, 

12 The notation /i points out that this term is stable, i.e. the algorithm will not change it 
any more. 
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Now, consider the cubic term fa; we normalize it by solving C0Q12) = ""0/25 
and we get f% = fa ~ £0(^2)- If this is nonzero we would like to (attempt 
to) further eliminate some terms. We can do it by exploiting the higher order 
effects of a transformation with generator , and in order to avoid this affects 

fa 1 ' we have to require £ Ker(£o). From the general formula (14) we have 
then that with such a transformation 

fa - /f = fi X) + (16) 
We define the linear "higher homological operators" £& as 

Ck(a) :={/ fc ,a} (17) 

so that the above reads 

If ] = fP-Cx(hP) , (18) 



and we can eliminate from all terms which are in the range of M.\, defined 
as the restriction of C\ to Ker(£ )- Thus, if we choose the complementary 
spaces to the ranges considered as the orthogonal complement, we get fa £ 
[Ran(£o)]- L n [R&n(Mi)] x . 

To avoid any misunderstanding, we stress immediately that our procedure 
is based on these restricted operators Mk, see below, and not on the full higher 
homological operators £&. 

Higher order terms, fa and so on, will also change; again we will denote 
them, for the sake of notation, still as fa etc. not to introduce heavy notation. 

Let us, for the sake of clarity, go one step further and consider the quartic 
term fa. First we normalize it by a transformation generated by /13, obtaining 
a /g 1 "* £ Ker(£ f J"). Then we can still use a transformation generated by £ 

Ker(£o) to eliminate some terms without affecting fa; we obtain in this way a 

(2) . . (3) 

fa . But we can also still use a transformation generated by h\ provided this 

does not affect fa nor fa, i.e. provided h\ £ Ker(£ ) H Ker(£i); with this, we 
obtain 

/3 

i.e. we can eliminate also all the terms which are in the range of M.^-, the 
restriction of £2 to Ker(£ ) n Ker(£i). 

At this point the general procedure should be quite clear; in order to state 
precisely the result I should define some chains of spaces and operators^. 

• Let #(°) := V, and define ij(P +1 ) := Ker(£ ) n ... n Ker(£ p ). Notice that 



fa = /f - £2 K 3) ) , (19) 



13 Note that all of these spaces could be intersected w ith and operators restricted to, the 
Vfc, obtaining finite dimensional spaces and operators [p5[ . I will not do this to avoid a 
cumbersome notation, but in computer algebra implementation of this algorithm one would 
proceed in this way. 
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• Define the operators M p := C P \ H ( P ). 

• Let := V, and define := [Rshi(Mq)} 1 - n ... n [Ran(Mp)]- 1 . Notice 
that C FO). 

It will also be convenient to set a notation for relevant projection operators: 

• Denote by Xs the projection operator \s '■ V — > Ker(£ s ) 

• Denote by /i s the composition of projection operators /i s = Xs-i ' ••• ' Xo f° r 
s > 1; we set fi = J- 

• Denote by tt s the projection operator n s : V — > Ran(.M s ). 

Note that with these definitions, = /x p y, and M. p = C p o fi p . 

Definition. The dynamical system x = f(x) = X^fcLi fk(%) is m Poincare 
renormalized form (PRF) up to order N if G F( fc ^ for all k < N. In the 
formal limit iV — > c», we say the system is in PRF, tout court. 

Theorem. Any analytic or formal smooth dynamical system can be taken into 
PRF up to any desired order N by means of a sequence of Lie-Poincare trans- 
formations. 

(v) 

Proof. Proceeding as above, the transformation generated by n£_ p will take 
fi P) into 

fl p+1) = fi p) - M P (AW,) (20 ) 
If we choose as a (non-unique) solution to the higher homological equation 

Tp4 P) -M p (hf_ p ) = (21) 



i.e. as 



h%> = » p -m;-k p n p > (22) 



k—p 



e(p) 



where A4* is the pseudo-inverse to A4 P , then we arrive - after applying the full 
procedure to fx, ...,/at - to a system x — f*(x) — J^'kLo fk( x ) wnere fk( x ) = 
£ F (k) n Vk for all k < N _ A 

Remark. The proof is completely constructive, as (22) tells which generator 
should be chosen at each step (and this is of course the nontrivial content of the 
theorem). It can thus be transformed in an algorithm to be implemented on a 
computer. 

Remark. The algorithm can deal with systems having zero linear part, contrary 
to the standard normalization algorithm. For such an application in a symmetry 
context, see p7| . 

Remark. The algorithm is only based on the grading of resonant terms (reso- 
nant vector fields) by their homogeneity properties, and is therefore completely 
general. However, it fails to recognize the Lie algebraic properties of the set 
of resonant vector fields. Not surprisingly, taking these into account leads to a 
more powerful reduction, which we discuss below. © 
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Notice that this procedure, as discussed here, is completely formal. A 
discussion (certainly not final) of convergence of PRFs - and also of PRF for 
Hamiltonian systems, and of PRF in the presence of symmetry - is contained 



in |26| and in my contributions to the SPT workshops ||9j, 29 . A discussion of 
these matters is certainly incompatible with the limited scope of this paper. 

I would now like to briefly mention an example of application of PRFs to 
a concrete problem, i.e. the (formal) analysis at singular points of vector fields 
in the plane having purely rotational linear part; this case corresponds to the 
vector field being, in the first order approximation, an integrable hamiltonian 
vector field. A complete analysis of regular singular points of vector fields in 
the plane from the point of view of PRFs is contained in |36|, where analytic 
computations are described in detail, and completely explicit computations are 
pursued up to order six (order ten in a special case). 

Example. Let the dynamical system x — f(x) in R 2 have a singularity at the 
point Xq; we can always shift coordinates so that xq — 0. Let the linearization 
x = [(djf % )(xo)]xP = Ajxi of X at x Q = be given by 

A - (J -) (23, 
Then the standard normal form is given (in vector notation) by 

oo 

x = Ax + VJ r 2k (a fe J + b k A) x (24) 
fc=i 

where I is the identity matrix and r 2 = \x\ 2 = (x'f + x 2 ,). 

In general, ai,b\ are nonzero (but if the system is hamiltonian, all the a k 
vanish; this case is treated in ^?j); however, for the sake of generality, let 
us assume = for k < fi and = for k < v, with a M ^ and 6„ ^ 0. 

Then, following step by step the procedure described above, the PRF is 
given by 

= Aw + (r 2M a M , + r 4tJ -ct2fi) w for fx < v 

= (1 + r 2 ^/3 M ) Aw + (r 2 ^a M + r 4 ^a 2M ) w for fi = v (25) 
= (1 + r 2v /3 u ) \Aw + E^(r 2fe a fe )w for (i > v 

where otk,/3k are real numbers (in general different from the ak,bk)- 

This result was obtained in |3Jj ; however the computation (and the result 

reported) there is not correct in the case with fj, > v; the correct computation 

leading to (25) is given in Q . 

Notice that the standard NF is characterized by two infinite sequences of 

constants; the PRF is finite except in case /i > v, where it however depends on 

only one infinite sequence of constants. 
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5 Lie renormalization 



Let us look more closely at the set Q of nonlinear vector fields resonant with a 
given linear one, Xa = (Ax) l di. 

First of all, it is clear that they form a Lie algebra (in general infinite dimen- 
sional) under commutation. Actually this also belongs to a finitely generated 
module, as we now briefly discuss recalling a general characterization of vector 
fields in normal form |^(| |33], ^j, ^59] . 

We say that the differentiable function (p : R" — > R is an invariant for 
Xa if Xa((p) — 0. Denote by 1(A) the set of invariants for Xa which are 
meromorphic (that is, can be expressed as a quotient of analytic functions) in 
the x coordinates; and denote by I*(A) C 1(A) the set of analytic invariants 
for X A . 

Let G = C(A) be the centralizer of A in the algebra of n dimensional 
matrices; let its Lie algebra be spanned by matrices {K%, K,i} (we can always 
assume K\ — I, and take one of the K a to be A if this is nonzero; notice that 
d < n). We denote by X^ the vector fields corresponding to these, i.e. given 
in the x coordinates by X^ — (K a x) l di. 

Then it results [^6|, |3^, ^?], |6f| that the most general vector field W in Q 
can be written as 

d 

W = ^(x) X< a > (26) 

where n a (x) £ 1(A); this means that Q is contained in the module over 1(A) 
generated by G. 

Notice that the vector field W must be analytic in the x, while generic 
vector fields of the form ^i a (x)X^ are not such: Q is not the full G-generated 
module over 1(A) . 

In several cases it happens that Q has a more convenient structure, i.e. the 
fjL a in (26) can actually be taken to be in 2* (A), and not just in T(A). In this 
case we have Q — I*(A) <E> G, and we say that all the vector fields in Q are 
quasi-linear, or that we have a quasi-linear normal form. In particular, this 
is e.g. the case when A admits only one basic invariant. 

Call X a the algebra spanned by vectors which are written as X — s(x)X^ a ' 
with s € 1(A); call X* C X a the algebra spanned by vectors as above with 
s G I*(A) (this is the full module over I*(A) generated by X^). As seen 
before, in general we have X{ © ... © X^ C Q c X\ © ... © Xd, and in the 
quasi- linear case we actually have Q = X* © ... © X%. 

Consider now the commutation relations between elements ^i a (^i)X^ and 
ap('$)X(P' of the subalgebras X a and Xp respectively ; it is immediate to check 
that, with Vj := d/dlpi, we have 

- (o p (V^ a )X^(^))X^ + (n a a ) [XM,XW] . {27) 
Remark. Notice that when X^ = X A , by definition X^(tpi) = 0, and 
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[X^ a \ X^] — 0; thus the corresponding subalgebra Xp is always an abelian 
ideal in Q. 

Remark. Note also that, as obvious from the formula (27) above, the direct sum 
of subalgebras X ai @...@X aa is a subalgebra in Q if and only if {X^ ai \ X^ 01 "'} 
span a subalgebra in G. 

Assume now that Q is quasi-linear, and that we are able to determine a 
sequence of subalgebras T p C £/, each of them being the direct sum of X a 
subalgebras, such that J-q = Q and 

[ Q , F P ] = J>f i ; (28) 

if this terminates in zero we say that Q has a quasi-nilpotent structure. Notice 
that the factor algebras T p := T p / J- p +\ are in general not abelian. 

By the above remark, Q can have a (nontrivial) quasi-nilpotent structure 
only if G is nilpotent. The chain of subalgebras T p C Q can then be read 
off the descending central series G p of G; recall that the factor algebras j p = 
Gp/G p +i for this are abelian. The nonabelian subalgebras T p introduced above 
are therefore modules over 1(A) generated by abelian subalgebras j p of G. 

Assume now Q is quasi-nilpotent. In this case we can first work with gener- 
ators in Ti and simplify terms in Fi (e.g. by following the PRF algorithm within 
the set I\; this allows to work with more familiar projection and homological 
equations than if setting the problem in a completely Lie algebraic framework) , 
then consider generators in T 2 and simplify the corresponding terms being guar- 
anteed that Ti terms are not changed, and so on. 

We will call the reduced normal form obtained in this way the Lie renor- 
malized form (LRF); the Lie renormalization is thus a blend of Lie algebraic 
reduction along the chain of subalgebras identified by (28) and of PRF reduc- 
tion in each of the r Q subsets. 

Remark. In this case we are - roughly speaking - just using the nilpotent 
structure of (the finite dimensional group) G, rather than the one of (the infinite 
dimensional algebra) Q. 

Remark. Needless to say, this approach is particularly convenient when the T p 
are generated by a single element of G; it can however be applied also in other 
situations. 

Let us briefly discuss why we have considered the quasi-linear case rather 
than the general one. It should indeed be clear that the key feature is the 
assumption that we are able to identify T p satisfying (28), so that this kind 
of approach is not confined to the quasi-linear case. However, recall that we 
have to consider smooth vector fields as generators of our Lie-Poincare changes 
of coordinates, and that at each step we have to deal with smooth (partially 
renormalized) vector fields. As these are only a subset of X\ © ... © X& (and one 
identified by analytic, not geometric, conditions), checking these conditions at 
each step will in general be a prohibitive task. The problem is not present for 
quasi-linear normal forms, when we deal with the full module T*(A) G. 

Thus, in practice, the LRF approach is not always applicable; however, 
when it can be applied the computations required are particularly simple fl36|. 



13 



We also stress that the LRF (like the PRF) for a given system is in general not 
unique. 

We stress that the quasi-nilpotent structure is met in applications: e.g., 
it applies to any nontrivial two-dimensional case and several three-dimensional 
ones Q. More generally, it always (but not only) applies when there is only 
one basic invariant. We consider now again a simple example; for other explicit 
ones see p6[ . 

Example. Let us consider again the case considered at the end of the previous 
section. Proceeding as discussed in this section, we obtain the LRF 

w = Aw + (r 2 ^a f _ l +r 4fi a 2 ^)w + {^^r 2k (3}}j Aw 

where the sum should be considered as zero for fi < v. Notice that now we get 
a finite dimensional reduced normal form in all cases. 



6 Filtering of Lie algebras 

I will finally briefly recall (and reformulate, hopefully without distorting it to 
such an extent to make it completely different from the original version) the gen- 
eral approach employing the Lie algebraic structure of Q (recall this included 
nonlinear fields only) proposed by Broer and Takens [[l7| [ll| , and further devel- 
oped by other authors. 

The Broer procedure is always applicable, and is able produce a unique 
reduced normal form ^ ^ . Unfortunately, implementing it in practice seems to 
be quite hard, and this approach has been applied only to systems of very low 
dimension. On the other side, the ideas put forward by Broer are at the basis 
of many approaches - including the one depicted in previous sections - to the 
problem of further reduction of normal forms. 

One can consider in general terms the descending central series of Q jl9], |54| . 
We recall this is defined by Go = Q ancf^l Gk+i = [G, Gk]] as well known the 
factor algebras = Gk/Qk+i are abelian. It is then possible to eliminate terms 
by inner automorphisms of Q (that is, by acting on Q with vector fields in Q) 
proceeding along Gk, i-c. filtering the Lie algebra Q. This is relevant to our 
problem since Lie-Poincare transformations always reduce to inner operations 
on Q. 

This approach is, of course, particularly convenient when Q is nilpotent; 
actually if Gk is the descending central series for G, we have Gk H V m = for 
m < k. We can reduce to the more familiar case of algebras nilpotent of finite 
order N by considering only the projections of Gk and Tk to the linear space 
V N = Vi © ... ®V N CV. 

Let us now briefly describe the implementation of the Broer renormaliza- 
tion procedure in the language employed in the present paper. We consider 

14 Needless to say, is not the set of fields in Q homogeneous of some degree in the a;! 
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the sequence of abelian factor algebras Tk = Gk/Gk+i, an d denote by Xk the 
projection from Q to IV 

We will produce a sequence of vector fields each being the result of 

the first j steps of further normalization. Let rjj be the operator on Q defined by 
rij(H) = e H W*^ e~ H ; let Wj be the operator of projection from Q to the range 
of rjj. We can reduce the normal form W = as follows. 

As the first step, consider a G Q and require that W = w!; 1) :=%(#) 
is such that noixoVoiH^)] = 0. This determines (non uniquely) H^°\ and 
produces a wi x \ 

At further steps, we have the same setting; the "homological equations" on 
Lie algebras to be solved at each step will be 

*i [Xj (vj(H^))] = 0; (29) 

this determines H^h Each wi J ' +1 ^ is then determined as 

Wi 3+1) = r)j{H {j) ) := exp[HW] exp[T^i j) ] cxp[-H^} . (30) 
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